QED theory of the nuclear magnetic shielding in hydrogen- like ions 
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The shielding of the nuclear magnetic moment by the bound electron in hydrogen-like ions is 
calculated ab initio with inclusion of relativistic, nuclear, and quantum electrodynamics (QED) 
effects. The QED correction is evaluated to all orders in the nuclear binding strength parameter 
and, independently, to the first order in the expansion in this parameter. The results obtained lay the 
basis for the high-precision determination of nuclear magnetic dipole moments from measurements 
of the ^-factor of hydrogen-like ions. 
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Magnetic dipole moments of nuclei are most often de- 
termined by the nuclear magnetic resonance (NMR) tech- 
nique. Other methods such as atomic beam magnetic res- 
onance, collinear laser spectroscopy, and optical pumping 
(OP) have also been used. The measured quantities are 
usually the ratio of the frequencies (or the ^-factors) for 
the nucleus of interest and the reference nucleus. Such 
ratios can be experimentally determined with a part-per- 
billion (ppb) accuracy However, magnetic moments 
of bare nuclei extracted from these experiments are much 
less accurate. This is because the experimental data 
should be corrected for several physical effects, which are 
difficult to calculate. The main effect is the diamagnetic 
shielding of the external magnetic field by the electrons 
in the atom. The NMR results should be also corrected 
for the paramagnetic chemical shift caused by the chem- 
ical environment Q and the OP data are sensitive to the 
hyperfine mixing of the energy levels [3| . Significant (and 
generally unknown) uncertainties of calculations of these 
effects often lead to ambiguities in the published values 
of nuclear magnetic moments Q. 

Since the accuracy of calculations of the chemical shifts 
cannot be reliably assessed, the means of comparison of 
nuclear moments shielded by different environments in 
NMR measurements are rather limited. Independent de- 
terminations of nuclear magnetic moments would define 
uncertainties of theoretical calculations of the chemical 
shifts and help to assess the accuracy of NMR standards. 

Reliable determination of the nuclear magnetic mo- 
ments is also prompted by a new generation of QED 
calculations of the hyperfine splitting in highly charged 
ions. It was demonstrated |5j that the magnetic sector of 
bound-state QED can be tested in these systems to all 
orders in the binding field, if the nuclear magnetic mo- 
ments are accurately known. Alternatively, comparing 
theoretical predictions with experimental results, one can 
determine nuclear properties and set benchmark tests for 
nuclear-structure theory. A recent example is the spec- 
troscopic determination of the nuclear charge radii of the 
neutron-halo nuclei 8 He, n Li, and n Be [6], which yielded 
unique information about the properties of these extraor- 
dinary systems. 



A way to a high-precision determination of nuclear 
magnetic moments is to study the simplest atomic sys- 
tems, the hydrogen-like ions. Measurements of the 
bound-electron ^-factor in these systems progressed dra- 
matically during the recent years and reached the ppb 
level [7]. They led not only to a stringent test of sophis- 
ticated QED calculations [H, ^ but also to an improved 
determination of the electron mass [lOj. Extensions of 
these experiments to ions with a nonzero nuclear spin 
will provide a determination of the nuclear magnetic mo- 
ments from a simple system that can be described theo- 
retically up to a very high accuracy. 

It is well known that the nuclear-spin-dependent 
part of the atomic ^-factor gp is suppressed by about 3 
orders of magnitude as compared to the leading effect 
due to the bound-electron ^-factor (see Eq. (pQ) below). 
This imposes limitations on possible determinations of 
the nuclear magnetic moment from gp. We show here, 
however, that the leading effect cancels exactly in the 
sum of the ^-factors for two hyperfine-structure levels 
(see Eq. (j2j) below). This sum is proportional to the 
nuclear ^-factor and, therefore, is much better suited for 
extracting the nuclear magnetic moment. Its calculation 
can be conveniently parameterized in terms of the nuclear 
shielding constant a, as given by Eq. (|2]). 

In this work we perform an ab initio calculation of 
the nuclear magnetic shielding for the ground state of 
hydrogen-like ions. The relativistic, QED, and nuclear 
effects are accounted for. The main challenge is the cal- 
culation of the QED correction. To the best of our knowl- 
edge, the only attempt to address it was the estimate re- 
ported in Ref. [l2j |. In this Letter, we calculate the QED 
correction rigorously to all orders in the binding nuclear 
strength parameter Za (where Z is the nuclear charge 
and a is the fine- structure constant) and, independently, 
we derive the leading term of its Za expansion. 

We now turn to the theory of the ^-factor of a 
hydrogen-like ion with a nonzero spin. Within relativistic 
quantum mechanics, it is given by 
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FIG. 1: Self-energy correction to the nuclear magnetic shielding. Double line represents the electron in the binding nuclear 
field. Wave line terminated by a triangle represents the dipole hyperfine interaction with the nucleus and wave line terminated 
by a cross represents the interaction with the external magnetic field. 



where F is the total angular momentum, / is the nu- 
clear spin, j is the electron angular momentum, gj is 
the Dirac bound-electron ^-factor, gj = /i/(/ijv/) is 
the nuclear ^-factor, /i is the nuclear magnetic mo- 
ment, fi 7v = \ e \/{^ m p) is the nuclear magneton, m and 
m p are the electron and proton masses, respectively, 
(j • F) = [F(F + 1) - 1(1 + 1) + j(j + l)]/2, and 
(I ■ F) = [F(F + 1) + 1(1 + 1) - j(j + l)]/2 . The 
higher-order corrections enter into Eq. (pQ) in two ways: 

(i) the Dirac electron ^-factor gj is modified by QED 
and recoil effects that do not depend on nuclear spin, 

(ii) the free-nucleus ^-factor gj is shielded by the bound 
electron. Additional corrections, e.g., those due to the 
electric quadrupole interaction [13j , are small and can be 
absorbed into the definition of the nuclear shielding. 

For the ground state of an ion with a nuclear spin / > 
1/2, we introduce the combination of ^-factors g, 

9 = 9f=i+i/2 + 9 F =1-^/2 = -2 (1 - a) , (2) 

m p fj, N I 

whith (7 being the shielding constant. If both (/-factors 
gF=i+i/2 an d 9f=i-i/2 are measured and a is known from 
theory, the above formula determines the nuclear mag- 
netic moment \i. For the ions with a nuclear spin / = x /2, 
Eq. (j2j) is not applicable and the nuclear magnetic mo- 
ment has to be determined from Eq. (pQ). 

The nuclear shielding constant a defined by Eq. (J2j) 
can be represented as a sum 

a = a (0) + SCTQED + S(T rec + S(J BW + 5(Tq , (3) 

where the leading-order relativistic result (includ- 

ing the finite nuclear size effect), ^ctqed is the QED cor- 
rection, 5<j rec is the recoil correction, Sctbw is the nuclear 
magnetization distribution (Bohr-Weisskopf) correction, 
and 5cfq is the electric quadrupole correction. 

The exact relativistic result for the leading-order mag- 
netic shielding was obtained analytically (for a point 
nucleus) 11411 and numerically [13] . The recoil correction 
is known [12| to the leading order in Za, 
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where M is the nuclear mass and gw = Mfi/ (hn I Z m p ). 
The exact relativistic result for the electric-quadrupole 
correction is [13| 
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where Q is the n uclear electric quadrupole moment and 

7 = \A-(W. 

We now turn to the QED correction to the nuclear 
magnetic shielding. It consists of the self-energy (SE) 
and vacuum-polarization parts, the SE being the most 
difficult one. The Feynman diagrams representing the 
SE correction (Fig. [T]) contain two magnetic interactions, 
one with the external magnetic field (in what follows, 
the Zeeman interaction), V zee (r) = ^ B ■ (fx a), and 
the other with the magnetic dipole nuclear field (in what 
follows, the hfs interaction), Vhf s (r) = ^ /x • (fx a)/r 3 . 
Formal expressions for the corresponding energy shifts 
can be obtained by the two-time Green's function method 
[l5j]. Irreducible parts of the diagrams in Fig. DJa)-(c) 
give rise to the perturbed orbital contribution, 

SE po = 2 (a\^(e a )\S^a) + 2 (S^a\Jl(e a )\5^a) , (6) 

where E is the SE operator, \silla) and |^fs a ) are the 
first-order perturbations of the reference-state wave func- 
tion induced by V zee and Vhf s? respectively, and \S^a) is 
the second-order perturbation induced by both interac- 
tions. The SE operator is defined by 
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where I(uj) = e 2 a^a u D^ u D^ v (uj) is the photon 
propagator, and u = 1 — iO. The diagram in Fig. QJd) 
gives rise to the hfs-vertex contribution, 

SE vrM s = 2 (a|r h fs(ea)|^> + 2 (a|S 7 |^a) (T4fs) , (8) 

where the prime denotes the derivative of the operator 
with respect to the energy argument and 
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The diagram in Fig. QJe) induces the Zeeman- vert ex con- 
tribution, in analogy with its hfs-vertex counterpart, 

5E vr , zee = 2 (a\T zee \5l)la) + 2 (a\E f \5^a) (V zee ). (10) 

Finally, Fig. [TJf ) together with the remaining derivative 
terms yields the double-vertex contribution, 

SE dmVt = 2(A) + (S // )(Kee)(^hfs) + (r[ lfs )(Kee) 

HKee)(VMs) + 2 (E 7 ) (a\V zee \5^a), (11) 
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TABLE I: QED corrections to the nuclear magnetic shielding. 
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SE 


VP 


10 


-0.51 (10) 


0.229 
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-0.710(15) 


0.256 
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-0.789 (9) 


0.271 
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-0.927 (4) 


0.302 
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-1.110(2) 
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-1.283 (1) 


0.417 


40 


-1.519(1) 


0.520 


54 


-2.029 (1) 


0.775 
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-4.457 (2) 


1.996 
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-7.107(2) 


2.954 
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FIG. 2: (Color online) Individual contributions to the nuclear 
shielding. "NR" is the nonrelativistic contribution, "REL" 
is the relativistic point-nucleus contribution, "FNS" is the 
finite nuclear size correction, "QED" is the QED correction, 
"BW" is the Bohr-Weisskopf correction, "REC" is the recoil 
correction, and "QUAD" is the electric quadrupole correction. 
Note that the QED correction changes its sign between Z — 4 
and 5. 



where A = A(e a ) is the 4-point vertex operator, 
(i\A(e)\k) =— dco J2 
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The formulas reported so far refer to the energy shifts. 
The corrections to the magnetic shielding are related to 
the energy shifts by Sa { = SEi IF(F + M F (I • 

F)), where Mp is the projection of the total momentum 
F. It can be shown that for the j = x j2 reference states, 
£ctqed does not depend on nuclear quantum numbers. 
The numerical calculation of c)<tse was performed along 
the lines developed in Ref. [16]; its details will be reported 
elsewhere. 

The remaining part of the QED effect is the vacuum 
polarization (VP). In our calculation, we include two 



dominant VP corrections induced by (i) modification of 
the electron line by the Uehling potential and (ii) modi- 
fication of the hfs interaction by the free-loop VP. 

Our calculational results for the SE and VP corrections 
are listed in Table U expressed in terms of the function 
D(Za), 



Sa Q EB = a 2 (ZafD(Za). 



(13) 



The SE correction is calculated for the point nucleus, 
whereas the VP part accounts for the finite nuclear size 
as well as higher-order iterations of the Uehling potential. 
Because of large numerical cancellations, we were able 
to perform our numerical SE calculations for Z > 10 
only. In order to extend our calculations to the lower- Z 
ions and to cross-check the numerical procedure, we also 
performed an analytical calculation of the leading term 
of the Za expansion. The result valid for an ns state 
reads 
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where In k (ls) = 2.984 128 and In k 3 (ls) = 3.272 806 @. 
Details of the analytical calculation will be reported else- 
where. The results of the numerical and the analytical 
calculations are in good agreement. 

We now turn to the effect induced by the spatial dis- 
tribution of the nuclear magnetic moment, also known as 
the Bohr-Weisskopf (BW) correction. Following Ref. [ItJ, 
our treatment of the BW effect is based on the effective 
single-particle model of the nuclear magnetic moment. 
Within this model, the magnetic moment is assumed 
to be induced by the odd nucleon with an effective g- 
factor, which is fitted to yield the experimental value of 
the nuclear magnetic moment. Under these assumptions, 
the BW effect can be described by the magnetization- 
distribution function F(r) that multiplies the standard 
point-dipole hfs interaction Vhfs( r )- The function F(r) is 
induced by the wave function of the odd nucleon, which 
is obtained by solving the Schrodinger equation with the 
Woods-Saxon potential (see Ref. [18] for details). The 
BW correction ^<jbw is obtained by reevaluating the 
leading-order magnetic shielding <j^ ^ with the hfs inter- 
action Vhf s multiplied by the magnetization-distribution 
function F(r). The relative uncertainty of 30% is as- 
cribed to this correction, which is consistent with previ- 
ous error estimates for this effect [l7| . 

Numerical results of our calculations are presented in 
Table HI and Fig. [2 The error of the QED correction 
comes from the numerical uncertainty of the SE part 
and the estimate of uncalculated VP terms (30% of the 
total VP part). The error of the quadrupole contribu- 
tion comes from the nuclear quadrupole moments. The 
largest error is due to the BW correction. Since this ef- 
fect cannot be presently accurately calculated, this uncer- 
tainty sets the practical limit to which the nuclear mag- 
netic moment can be determined from an atomic system. 
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TABLE II: Individual contributions to the shielding constant a x 10 6 for selected hydrogen- like ions, see Eq. ©• 



i7 7+ 43 Ca 19+ 73 Ge 31+ 131 Xe 53+ 209 Bi 82+ 



Leading 143.3127 375.960 657.93 1461.6 4112 

QED -0.0026(2) -0.103(15) -0.59(8) -4.1 (0.8) -30(7) 

Bohr-Weisskopf -0.0013(4) -0.061 (18) -0.54(16) -8.2(2.5) -42(13) 

Quadrupole -0.0007(1) -0.018 -0.42 6.9(0.1) 7 

Recoil -0.0120 -0.015 -0.02 0.0 

Total 143.2960(5) 375.763(24) 656.36(18) 1456.3(2.6) 4046 (15) 



For very light ions, the theoretical accuracy is limited 
by the recoil effect (see Fig. [2j), which is known in the 
nonrelativistic limit only. Note that some of corrections 
to a depend on the nuclear ^-factor. This dependence, 
however, is so weak that it can be safely ignored in the 
determination of the magnetic moments. 

Summarising, we have presented ab initio calculations 
of the nuclear shielding in hydrogen-like ions, which ac- 
count for relativistic, nuclear, and QED effects. The 
present theory permits determination of nuclear mag- 
netic moments with fractional accuracy ranging from 
10- 9 in the case of 17 7+ to 10" 5 for 209 Bi 82 +. This 
Letter is primarily focused on nuclei with spin / > x /2, 
but the case of / = x /2 is only slightly more complicated. 
Then, the nuclear-spin-independent part of gp in Eq. (pQ) 
can be cancelled approximately, by taking a difference of 
the factors gp for two different isotopes of the same 



element. 

Modern experiments on ^-factors of hydrogen-like ions 
have achieved the accuracy of a few parts in 10 11 [19[ but 
so far have been restricted to ions with spinless nuclei. 
Their extention to the nuclei with spin requires driving 
the hfs transition and measuring the ^-factor of an atom 
in a hyperfine excited state. These are significant compli- 
cations but they do not make an experiment prohibitively 
difficult 0. 
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